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Conversion  Factors,  Non-SI  to  Si 
Units  of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI  (System 
Internationale)  units  by  applying  the  following  factors: 


Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

miles  (U.S.  statute) 

1.6094 

kilometers 

tons 

907.185 

kilograms 

IX 


1  Introduction  and  Background 


This  report  describes  the  prediction  of  loading  cycles  on  miter  gates  for  use  in 
the  assessment  of  fatigue  reliability  for  miter  gates.  The  report  explains  the 
correlation  of  field  data  for  pool  and  tailwater  elevations  and  barge  traffic  to  form 
a  loading  histogram  to  be  utilized  to  better  predict  the  loading  history  of  a  miter 
gate. 

Miter  gates  at  navigation  locks  experience  loading  cycles  from  emptying  and 
filling  of  a  lock’s  chamber  as  they  are  opened  to  allow  traffic  through  the 
navigation  locks.  Reliability  analysis  of  miter  gates  at  navigation  locks  requires  the 
definition  of:  (a)  nonperformance  modes,  (b)  loads,  (c)  structural  strength,  and  (d) 
methods  of  reliability  analysis.  Due  to  the  cyclic-loading  nature  of  miter  gates,  the 
fatigue  of  critical  details  requires  examination  using  reliability  methods.  The 
assessment  of  fatigue  reliability  of  these  details  as  a  function  of  time  requires  the 
knowledge  of  strength,  stress  ranges,  and  loading  cycles  for  these  details. 

The  strength  of  fatigue  details  can  be  expressed  in  the  form  of  Stress-range 
versus  Number  of  cycles  (S-N)  curves.  The  determination  of  stress  ranges  requires 
analyzing  the  reliability  of  miter  gates  under  different  loading  conditions.  The 
number  of  loading  cycles  needs  to  be  determined  also  because  it  constitutes  an 
important  component  in  the  reliability  analysis,  in  addition  to  its  use  in  determining 
the  reliability  as  a  function  of  time.  As  a  result  of  the  analysis,  an  evaluation  of  the 
remaining  life  of  critical  fatigue  details  in  miter  gates  can  be  made.  The  accuracy 
of  the  assessed  remaining  life  depends  in  part  on  the  accuracy  in  the  utilized 
loading  cycles  in  the  reliability  analysis. 


Objectives 

The  study  that  is  reported  herein  had  the  following  objectives: 

a.  Reviewing  previous  U.S.  Army  Corps  of  Engineers  (USAGE)  studies  for 
determining  loading  cycles  for  miter  gates  at  navigation  locks. 

b.  Determining  the  factors  that  affect  the  loading  cycles  of  miter  gates  at 
navigation  locks. 
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c.  Developing  a  methodology  for  computing  the  number  of  loading  cycles  of 
miter  gates  at  navigation  locks. 

d.  Assessing  future  traffic  using  the  General  Equilibrium  Model  (GEM)^  and 
relating  its  results  to  loading-cycles  number  for  miter  gates. 

e.  Demonstrating  the  use  of  the  GEM  for  assessing  future  traffic,  and  in 
relating  its  results  to  the  number  of  loading  cycles. 

/.  Demonstrating  the  methodology  using  a  case  study. 


Literature  Review 


The  fatigue  failure  mode  is  commonly  considered  in  detail  design,  after 
principal  structural  members  have  been  sized.  Several  procedures  have  been  used 
for  the  assessment  of  fatigue  damage  (Wirsching  1984,  Wirsching  and  Chen  1987), 
such  as  deterministic  methods.  Spectral  method,  WeibuU  model,  and  Nolte- 
Hasford  model.  For  example,  the  Spectral  method  was  used  for  marine  structures. 
As  was  demonstrated  by  Chen  and  Mavrakis  (1988),  the  Spectral  method  is  more 
accurate  than  the  WeibuU  model  for  the  case  of  offshore  platforms  because  its 
results  are  less  sensitive  with  respect  to  the  variability  in  the  shape  of  the  wave 
spectra  compared  to  the  results  of  the  WeibuU  model.  However,  the  Spectral 
method  is  also  the  most  computationaUy  intensive.  Fatigue  reUabUity  can  be 
evaluated  by  using  Munse’s  model  (Munse  et  al.  1982),  Wirsching’s  model 
(Wirsching  1984),  or  advanced  second  moment  methods  (Madsen,  Skjong,  and 
Moghtaderi-Zadeh  1986).  A  reliabUity-based  design  format  for  fatigue  was 
demonstrated  by  White  and  Ayyub  (1987). 

Fatigue  can  result  from  the  cycUc  loading  appUed  to  miter  gates  of  navigation 
locks  as  they  are  operated  to  perform  lockages  of  vessels  going  through  the  locks. 
CycUc  loading  on  a  gate  is  generated  by  the  cycUc  head  and  taUwater  pressures  on 
the  gate.  Therefore,  variable  ampUtude  stress  ranges  are  appUed  to  critical  fatigue 
detaUs  of  the  gate.  These  stress  ranges  cause  damage  to  the  detaUs  that  is 
cumulative  in  nature  untU  a  crack  initiates,  then  propagates  to  a  possible  faUure  of 
the  affected  structural  members  of  the  gate.  Miner’s  rule  of  cumulative  fatigue 
damage  under  variable  ampUtude  stress  range  with  stress-range  versus  number  of 
cycles  to  faUure  (S-N)  information  about  fatigue  detaUs  can  be  used  to  assess  the 
reUabUity  of  the  details  (American  Society  of  CivU  Engineers  (ASCE),  Committee 
on  Fatigue  and  Fracture  ReUabUity  1982;  White  and  Ayyub  1987;  Rides  and  Leger 


’  For  convenience,  symbols  and  abbreviations  are  listed  in  the  notation  (Appendix  B). 
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1993;  Sommer,  Nowak,  and  Thoft-Christensen  1993;  USAEWES  1994).  The 
details  of  fatigue  analysis  of  miter  gates  are  provided  in  ETL  1 1 10-2-346 
(Headquarters,  Department  of  the  Army  1993a).  Also,  ETL  1 1 10-2-351 
(Headquarters,  Department  of  the  Army  1993b)  covers  fatigue  analysis  of  spillway 
gates. 

Fatigue  reliability  requires  the  definition  of  the  S-N  curves  for  the  fatigue 
details  with  their  associated  uncertainties,  and  the  cyclic  variable  amplitude  stress 
ranges  with  their  modeling  uncertainty  in  the  form  of  stress  range  histograms 
(Ayyub  and  White  1990;  Ayyub,  White,  and  Purcell  1989;  Ayyub  et  al.  1990; 
USAGE  1994;  White  and  Ayyub  1987).  Therefore,  the  data  needed  for  fatigue 
reliability  assessment  are  (a)  types  of  fatigue  details,  (b)  statistical  description  of 
the  S-N  curves  for  the  fatigue  details  (Fisher  et  al.  1970, 1974),  (c)  materials  and 
their  strength  and  stiffiiess  properties  such  as  the  modulus  of  elasticity  and  Poisson 
ratio,  (d)  structural  geometry  and  dimensions,  (e)  joint  histograms  of  loading 
cycles  and  water  heads  on  both  sides  of  the  gate  as  a  function  of  time,  and  (f) 
modeling  uncertainties  in  Miner’s  rule  and  computed  stresses.  The  applied  loads 
can  also  include  impact  and  vibration  loads.  Their  effect  on  fatigue  depends  on  the 
occurrence  frequency  of  their  stress  ranges.  Several  recent  USAGE  studies  on 
fatigue  and  concrete  deterioration  include  information  on  loading  cycles  of  miter 
gates  at  navigation  locks.  These  studies  were  examined  for  any  information  and 
techniques  that  pertain  to  this  study  (Headquarters,  Department  of  the  Army 
1993a;  Headquarters,  Department  of  the  Army  1993b;  USAEWES  1994). 
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2  Factors  Affecting  Loading 
Cycies 


The  number  of  loading  cycles  for  miter  gates  is  a  random  variable  with  inherent 
uncertainty  resulting  from  sources  that  include  the  following: 

a.  Navigation  traffic  volume. 

b.  Traffic  composition  primarily  in  terms  of  tow  length. 

c.  Length  and  capacity  of  navigation  locks. 

d.  Traffic  direction  and  pattern. 

e.  Weather-related  conditions,  e.g.,  ice  buildup  and  debris. 

/.  Impact  loads. 


The  loading  cycles  need  to  be  expressed  in  a  form  that  is  suitable  for  computing 
variable-amplitude  stress  ranges  with  associated  loading  cycles.  Therefore,  the 
loading  cycles  should  be  related  to  pool  and  tailwater  elevations. 

The  navigation  traffic  volume  is  an  important  factor  in  determining  the  number 
of  cycles.  The  number  of  cycles  is  increased  by  increasing  the  traffic  volume. 
However,  the  relationship  is  complicated  by  the  traffic  composition  in  terms  of 
tow-length  distribution.  Tows  are  typically  longer  than  most  navigational  locks, 
and  are  passed  through  a  lock  in  several  sections  or  cuts  and  reassembled 
afterwards.  The  number  and  length  of  the  cuts  in  a  tow  is  dependent  upon  the 
lock’s  dimensions. 

Traffic  in  navigation  locks  can  be  in  either  an  upstream  or  downstream 
direction  or  both.  In  cases  where  the  traffic  includes  the  simultaneous  upstream 
and  downstream  movements,  the  gates  can  be  operated  more  efficiently  by 
alternating  between  the  two  types  of  traffic,  hence  utilizing  each  loading  cycle  of 
the  gates  in  moving  the  traffic. 

Winter  weather  conditions  affect  the  loading  cycles,  especially  ice  buildup  in  the 
upper  approach  regions  of  locks.  Some  locks  are  completely  closed  to  river  traffic 
during  the  winter  because  of  ice  buildup  on  the  river;  other  locks  are  operated 
year-round.  Sometimes,  miter  gates  are  operated  for  the  purpose  of  passing  ice 
flows  to  reduce  ice  buildup  in  the  upper  approach,  and  to  relieve  any  pressure  on 
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the  gates.  The  loading  cycles  for  managing  the  ice  flow  are  usually  not  recorded  in 
operational  logs  (Patev  1995). 

Sometimes  barges  or  boats  apply  a  load  to  miter  gates  due  to  unintentional 
impact  resulting  from  poor  judgment  of  tow  operators,  poor  hydraulic  conditions, 
or  inclement  weather.  Even  though  these  impact  loads  can  cause  damage  to  the 
gates,  they  are  not  considered  in  this  report  for  assessing  fatigue  reliability. 
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3  Information  Sources  for 
Assessing  Loading  Cycles 


Lock  Performance  Monitoring  System  (LPMS) 

The  USAGE  Lock  Performance  Monitoring  System  (LPMS)  includes 
information  such  as  lock  number,  lockage  date,  start  of  lockage,  direction  of 
lockage,  number  of  cuts,  entry  type,  exit  type,  end  of  lockage,  and  tonnage. 

LPMS  entries  can  be  used  to  compute  loading  cycles  needed  for  fatigue-reliability 
evaluation.  However,  the  computation  of  these  cycles  from  the  LPMS  can  require 
a  significant  level  of  effort  due  to  the  structure  of  LPMS. 


Lockmaster  Records 

Lockmasters  maintain  records  that  might  contain  information  on  loading  cycles 
of  miter  gates.  For  this  study  several  navigation  locks  were  selected  for  a  field  trip 
after  consultation  with  USAGE.  A  field  trip  to  the  sites  of  these  locks  was 
undertaken.  The  records  of  the  locks  were  examined.  Information  that  can  be 
used  in  developing  the  methodology  for  this  study  was  gathered.  The  field  trip 
included  Locks  and  Dams  22,  24,  and  25  on  the  Mississippi  River. 


Hydraulic  Records 

Hydraulic-operation  records  include  pool  and  tailwater  elevations  that  can  be 
obtained  on  daily  basis.  These  records  are  needed  for  assessing  loading  cycles 
with  associated  water  heads.  The  USAGE  records  were  found  to  be  suitable  for 
water-head  load  estimation.  The  pool  elevation  (or  height)  of  water  {Hp)  and  the 
tail  elevation  (or  height)  of  water  (ET,)  are  needed.  These  quantities  need  to  be 
computed  on  a  daily  basis  starting  from  the  beginning  of  the  LPMS  (i.e.,  January 
1980)  to  the  present.  The  heights  (HpJI,)  can  be  used  to  compute  stresses  and 
stress  ranges  at  critical  fatigue  locations,  whereas  the  number  of  repetitions  of  the 
pairs  (Hp,Ht)  produces  the  needed  frequency  of  the  corresponding  stress  ranges. 
Therefore,  stress-range  frequency  histograms  necessary  for  fatigue  analysis  can  be 
produced.  The  number  of  repetitions  of  the  pairs  (Hpji,)  can  be  computed  from 
data  in  the  LPMS. 
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Uncertainty  in  Loading  Cycles 


The  number  of  loading  cycles  for  miter  gates  is  considered  in  this  study  to  be  a 
random  variable.  The  sources  of  uncertainty  in  assessing  loading  cycles  are  due  to 
the  factors  discussed  previously.  The  uncertainty  in  the  number  of  loading  cycles 
can  be  expressed  using  a  probability  distribution  for  this  number.  Also,  the  water 
elevations  with  the  associated  loading  cycles  need  to  be  expressed  in  probabilistic 
terms.  Prediction  models  that  can  be  developed  based  on  statistical  analyses  of 
data  include  statistical  variability  that  needs  to  be  assessed.  Modeling  uncertainty 
that  is  associated  with  computing  stresses  at  critical  fatigue  locations  of  miter  gates 
needs  also  to  be  considered  in  performing  fatigue  reliability  assessment. 
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4  Methodology  for  Determining 
Loading  Cycles 


Stages  for  Computing  Loading  Cycles 

The  methods  needed  for  determining  loading  cycles  for  miter  gates  is 
dependent  upon  the  nature  of  available  information.  The  type  of  available 
information  depends  on  the  years  (i.e.,  time  period)  of  interest.  Loading  cycle 
estimates  should  cover  the  entire  life  of  a  gate  since  these  estimates  are  needed  for 
fatigue  analysis.  The  life  of  the  gate  can  be  viewed  to  consist  of  its  present  age 
plus  the  planned  remaining  life.  Estimates  of  loading  cycles  up  to  the  present  can 
be  determined  using  methods  that  depend  on  the  type  of  available  information. 
Estimates  of  future  loading  cycles  need  to  be  based  on  forecasting  models  of  future 
traffic  in  navigation  channels. 

In  this  chapter,  the  terms  lockage,  lockage  cut,  hydrostatic  loading  cycles,  and 
hardware  cycles  are  used.  In  general,  a  lockage  is  defined  as  a  series  of  events 
required  to  transfer  a  tow  or  vessel  with  all  its  barges  through  a  lock  in  a  single 
direction.  For  the  purposes  of  this  report,  a  lockage  cut  is  defined  as  a  process  of 
passing  one  cut  of  a  tow  or  several  vessels  together  through  a  lock.  This  process 
requires  the  operation  of  the  gates  of  the  lock  (the  emptying  and  filling  of  a  lock’s 
chamber)  once,  if  the  gates  are  favorably  positioned  to  an  inbound  cut  of  a  vessel. 
If  a  vessel  can  be  accommodated  in  the  lock  in  its  entirety,  then  one  emptying  and 
one  filling  of  the  lock’s  chamber  are  required.  Also,  if  several  vessels  can  be 
simultaneously  accommodated  in  the  lock,  then  one  emptying  and  one  filling  of  the 
lock’s  chamber  are  required.  However,  if  a  vessel  is  too  large  to  be 
accommodated  in  the  lock,  it  is  separated  in  two  or  more  cuts.  Several  lockage 
cuts  in  this  case  are  required  in  order  to  pass  through  all  the  cuts  by  emptying  and 
filling  the  lock’s  chamber  several  times.  The  number  of  lockage  cuts  in  this  case  is 
equal  to  the  number  of  cuts  passed  in  the  lock.  If  a  lock’s  state  is  not  in  a 
favorable  position  to  receive  an  inbound  vessel,  an  additional  cycle  of  emptying 
and  filling  of  the  lock’s  chamber  is  required.  A  hydrostatic  loading  cycle  consists 
of  a  complete  emptying  and  filling  of  a  lock’s  chamber  that  produces  a  hydrostatic 
water-head  differential  on  the  gates.  A  hardware  cycle  is  a  complete  emptying  or 
filling  of  a  lock’s  chamber  that  produces  a  hydrostatic  water-head  differential  on 
the  gates.  Therefore,  a  hydrostatic  loading  cycle  consists  of  two  hardware  cycles. 
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A  key  source  of  information  for  estimating  loading  cycles  on  miter  gates  is  the 
USAGE  LPMS.  The  LPMS  generally  covers  the  period  from  1980  to  present. 

The  time  period  before  1980  can  be  broken  down  into  several  stages  depending  on 
the  available  information.  For  example,  it  can  be  broken  down  into  two  stages, 
start  of  life  (completion  of  construction,  e.g.,  around  late  1930’s  and  early  1940’s 
for  most  locks  on  the  Mississippi  River)  to  1948,  and  1948  to  1980.  It  seems  that 
formal  traffic  record  keeping  for  these  locks  was  not  established  until  the  late 
1940’s. 

In  consideration  of  the  above  description  of  the  nature  of  available  information 
on  the  utilization  of  locks,  the  following  stages  can  be  identified  for  developing 
methods  for  estimating  loading  cycles: 


Stage  1:  Start  of  life  (e.g.,  1940)  to  1948. 

Stage  2:  1948  to  start  of  the  LPMS  (i.e.,  1980). 

Stage  3:  Start  of  LPMS  (1980)  to  present. 

Stage  4:  Present  to  planned  design  (or  rehabilitation)  life. 


In  this  chapter,  methods  for  assessing  loading  cycles  for  these  stages  were 
developed.  Stage  3  (i.e.,  the  LPMS  stage)  contains  information  of  the  highest 
levels  of  data  quality  and  certainty.  Therefore,  this  stage  can  be  used  as  a  basis  for 
estimating  some  of  the  parameters  in  other  stages  as  described  below. 


Loading  cycles  for  Stage  3  (start  of  LPMS  to  present) 

This  stage  contains  the  best  and  needed  information  to  evaluate  loading  cycles 
on  miter  gates.  The  primary  source  of  information  is  the  LPMS;  therefore,  this 
stage  is  called  hereafter  the  LPMS  stage. 

In  the  LPMS  stage,  the  following  quantities  are  of  interest:  (a)  pool  elevation 
(or  height)  of  water  {Hp),  (b)  tail  elevation  (or  height)  of  water  (//,),  and  (c)  the 
corresponding  number  of  repetitions  of  the  pair  (HpJh).  Due  to  the  observed  daily 
variability  in  the  water  elevations,  these  quantities  need  to  be  computed  on  a  daily 
basis  starting  from  January  1980  to  present.  The  heights  (HpJHl,)  can  be  used  to 
compute  stresses  and  stress  ranges  at  critical  fatigue  locations,  whereas  number  of 
repetitions  of  the  pairs  iHp,H,)  produces  the  needed  frequency  of  the 
corresponding  stress  range.  Therefore,  stress-range  frequency  histograms  that  are 
necessary  for  fatigue  analysis  are  produced.  The  daily  upper  and  lower  water  pool 
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heights  can  be  obtained  from  the  hydraulic  records  of  a  lock  as  discussed  in 
Chapter  3. 

The  number  of  repetitions  of  the  pairs  can  be  computed  from  data  in 

the  LPMS.  The  fields  of  the  LPMS  that  are  shown  in  Table  1  can  be  utilized  for 
this  purpose.  These  fields  are  defined  in  the  LPMS  user’s  manual  (USAGE  1990). 
The  number  of  hardware  cycles  in  a  day  (LR-SHFT-DY)  for  a  selected  month 
(LR-SHFT-MO),  a  selected  year  (LR-SHFT-YR),  and  a  selected  lock  (LR-LOCK) 
can  be  computed  from  Table  1  by  considering  for  the  start  of  lockage,  end  of 
lockage,  entry  type,  exit  type,  vessel  type,  and  direction  of  traffic  in  its 
computation.  Several  algorithms  were  developed  to  compute  hardware  cycles 
from  the  fields  in  Table  1.  Some  of  these  algorithms  produced  erroneous  results 
due  to  some  ambiguity  in  the  meaning  of  some  fields  of  the  LPMS.  Considering 
these  trials,  some  fields  of  the  LPMS  can  be  improved  to  facilitate  the 
computations  of  hardware  cycles.  For  example,  the  following  observations  are 
made: 

a.  The  current  entry  and  exit  types  do  not  necessarily  reflect  the  turnback 
type,  if  it  was  delayed,  i.e.,  not  immediate  to  an  entry  or  exit,  respectively. 
Depending  on  the  use  of  these  fields  in  their  current  forms,  either  new 
fields  should  be  developed  that  correct  for  the  delayed  turnback 
occurrence,  or  the  current  fields  should  be  revised. 

b.  Sometimes  several  vessel  records  were  entered  in  the  LPMS  as  separate 
lockages,  but  these  vessels  were  serviced  in  the  same  operation  of  opening 
and  closing  of  miter  gates.  The  LPMS  does  not  keep  track  of  these  cases, 
thereby  complicating  the  computation  of  hardware  cycles.  The  start  of  a 
lockage  and  end  of  a  lockage  were  used  to  account  and  correct  for  this 
factor. 

c.  Ice  and  debris  lockages  are  not  included  in  the  LPMS.  The  results  of 
time-lapsed  videotapes  of  Locks  22  and  25  were  used  to  assess  these 
cycles  (Patev  1995). 

d.  The  operations  of  the  gates  for  service,  inspection,  or  performance 
evaluations  are  not  recorded  in  the  LPMS. 


Considering  the  above,  the  daily  numbers  of  hardware  cycles  were  computed 
using  the  following  logic: 

a.  The  number  of  cuts  NC  (Le.,  LR-NO-CUTS  field)  was  corrected  to 
account  for  cases  that  involve  several  vessels  serviced  by  the  same 
lockage.  Therefore,  a  new  field  was  added  to  Table  1  called  the  corrected 
number  of  cuts  (CNC).  The  entry  of  this  field  for  the  fth  record  in  a 
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month  was  computed  using  the  following  logic  statement  in  spreadsheet 
form: 

CNC(0  =  IF(AND(DY(i)  =  DY(/+1) ,  OR(ABS(SOLl(f+l)  -  SOLl(O) 

<=  13,  ABS(EOLl(i+l)  -  EOLl(O)  <=  10) ,  NC(/+1)  =  NC(/) ,  VT(0  <> 
“R”) ,  NC(/)  - 1 ,  NC(0)  (1) 

The  variables  in  Equation  1  are  defined  in  Table  1,  where,  for  example, 
DY(0  =  day  of  shift  for  the  zfh  record  in  a  month,  and  VT(0  <>  “R” 
means  the  vessel  type  of  this  record  does  not  equal  “R”  which 
corresponds  to  a  recreational  vessel.  Equation  1  is  based  on  an  EXCEL 
(Microsoft®)  logical  statement  of  the  following  type: 

TF{logical  expression  ,  expression  (fTRUE  ,  expression  z/ FALSE)  (2) 

b.  Another  new  field  was  introduced  called  hardware  cycles  (HC)  for  the  ith 
record  in  a  month.  This  field  was  computed  using  the  following  logic 
statement: 

HC(f)  =  2*CNC(i)  - 1  +  IF(DR(0  =  DR(i+l)  ,1,0)  (3) 

c.  The  daily  number  of  hardware  cycles  was  then  computed  by  totaling  the 
values  of  HC(0  over  all  i  values  (Le.,  records)  that  are  in  the  same  day. 

d.  The  daily  number  of  lockages  (LG)  was  then  computed  by  totaling  the 
values  of  CNC(0  over  all  i  values  (i.e.,  records)  that  are  in  the  same  day. 

e.  The  daily  hardware-cycle  and  lockage  numbers  were  then  corrected  by 
adding  cycles  needed  to  pass  through  ice  and  debris.  These  additional 
cycles  can  be  estimated  based  on  time-lapsed  videotapes,  if  available. 

/.  The  monthly  number  of  hardware  cycles  was  then  computed  by  totaling 
the  values  of  items  c  and  e  over  each  month. 

g.  The  monthly  number  of  lockages  was  then  computed  by  totaling  the 
values  of  items  d  and  e  over  each  month. 
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Table  1 

Selected  Fields  fronn  the  LPMS 

Field  Number 

Field  Name 

Description 

1 

LR-LOCK 

Lock  Number  (LN) 

7 

LR-SHFT-MO 

Month  of  Shift  (MO) 

8 

LR-SHFT-DA 

Day  of  Shift  (DY) 

g 

LR-SHFT-YR 

Year  of  Shift  (YR) 

10 

LR-SOL-1 

Start  of  Lockage  Time  (24  hr)  1st  Cut 
(SOLI) 

12 

LR-DIR 

Direction  of  Lockage  (up  or  down)  (DR) 

13 

LR-NO-CUTS 

Number  of  Lockage  Cuts  (NC) 

14 

LR-LCKG-TYPE 

Lockage  Type  (LT) 

15 

LR-VSL-TYPE 

Vessel  Type  (VT) 

18 

LR-ENTRY-TYPE 

Entry  Type  (ET) 

19 

LR-EXIT-TYPE 

Exit  Type  (XT) 

26 

LR-EOL-1 

End  of  Lockage  Time  (24  hr)  First  Cut 
(EOL1) 

31 

LR-EOL-2 

End  of  Lockage  Time  (24  hr)  Last  Cut 
(EOL2) 

78 

LR-TONS 

Tonnage  (TN) 

The  logic  above  is  based  on  the  assumption  that  the  idle  gate  position  of  a  lock 
is  its  position  at  the  end  of  the  previous  lockage.  The  logic  above  accounts  for  the 
additional  hydrostatic  loading  cycles  needed  to  position  the  gates  in  a  favorable 
position  for  receiving  incoming  vessels  to  the  lock.  This  effect  was  accounted  for 
by  considering  the  sequences  for  the  direction  of  traffic.  The  entry  and  exit  types 
recorded  in  the  LPMS  were  found  to  be  unsuitable  for  this  purpose  because  these 
entries  show  only  the  immediate  entry  or  exit  types,  respectively. 

The  data  reduction  and  analysis  of  the  hydraulic  records  and  the  UPMS’s  fields 
produce  daily  quantities  for  pool  water  elevation  (Hp),  tailwater  elevation  (ff,), 
number  of  hardware  cycles  (HC),  and  number  of  lockage  cuts  (LGC).  These 
results  can  be  used  to  compute  the  loading  cycles  of  interest  by  following  the  steps 
below: 

Step  1.  Use  curve-fitting  to  develop  a  relationship  between  pool  water 
elevation  and  tailwater  elevation.  The  result  can  be  expressed  as  pool 
water  elevation  as  a  function  of  tailwater  elevation. 

Step  2.  Sum  the  numbers  of  hardware  cycles  for  intervals  of  tailwater 
elevations.  Normalizing  the  number  of  hardware  cycles  by  the  total 
number  of  hardware  cycles  produces  tailwater  elevations  with  associated 
fractions  of  hardware  cycles,  i.e.,  a  histogram  of  tailwater  elevation  in 
which  tailwater  elevation  is  treated  as  a  loading. 

Step  3.  Fit  a  probability  density  function  to  the  histogram  from  Step  2. 
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Step  4.  Determine  the  total  number  of  hardware  cycles  and  lockage  cuts  on 
monthly  and  yearly  bases  from  the  daily  records. 

Step  5.  Use  curve-fitting  to  establish  relationships  among  the  following 
variables:  hardware  cycles,  lockage  cuts,  tonnage,  and  time.  The  tonnage 
for  a  lock  over  some  time  period  can  be  computed  from  the  fields  of 
LPMS  as  the  total  weight  of  commodities  that  pass  through  the  lock 
within  this  time  period.  These  relationships  are  needed  in  other  stages  as 
discussed  below.  The  relationships  can  be  developed  using  monthly  or 
annual  records. 

The  models  that  result  from  Steps  3  and  5  constitute  the  basis  for  assessing  the 
loading  cycles.  These  models  can  be  expressed  in  dimensionless  format  by 
normalizing  them  with  respect  to  corresponding  design  values.  For  example,  the 
tailwater  loading  probability  density  function  (Step  4)  can  be  normalized  with 
respect  to  the  tailwater  elevation  design  value  for  an  investigated  lock.  The  benefit 
of  expressing  the  results  in  a  normalized  format  is  in  potentially  increasing  the 
range  of  applicability  of  the  results  to  other  locks.  In  general,  locks  and  dams 
along  a  river  can  be  classified  into  groups  (or  reaches).  A  typical  lock  can  be 
analyzed  from  each  reach  and  several  locks  in  a  selected  reach  can  be  analyzed  to 
produce  a  complete  understanding  of  loading  cycles  on  miter  gates.  Future  work 
in  this  area  can  examine  the  relationships  and  variability  among  the  reaches  and 
within  reaches. 


Loading  cycles  for  Stage  1  (start  of  life  to  1948) 

Stage  1  is  defined  as  the  stage  from  the  start  of  life  of  a  lock  to  about  1948. 
The  end  of  this  stage  (Le.,  1948)  was  established  based  on  the  record-keeping 
practices  of  the  USAGE  for  locks  on  the  Mississippi  River.  This  stage  can  be 
characterized  as  a  stage  with  inadequate  records  that  are  needed  to  compute 
loading  cycles.  Generally,  the  information  available  from  records  in  this  stage  is 
limited  to  annual  tonnage.  Therefore,  the  relationships  of  hardware  cycles  and 
lockage  cuts  as  functions  of  tonnage  that  were  developed  m  the  LPMS  (third) 
stage  can  be  used  to  estimate  hardware  cycles  and  lockage  cuts,  respectively,  for 
this  stage.  Alternatively,  the  number  of  hardware  cycles  as  a  function  of  time  and 
the  number  of  lockage  cuts  as  a  function  of  time  that  were  developed  in  the  LPMS 
stage  can  be  used  to  estimate  these  quantities  by  extrapolation.  The  former 
approach  is  recommended  in  this  report  because  its  results  are  partly  based  on  data 
in  the  first  stage. 
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Loading  cycles  for  Stage  2  (1948  to  start  of  LPMS) 

Stage  2  is  defined  as  the  stage  fi-om  about  1948  to  the  start  of  the  LPMS  stage. 
The  start  of  this  stage  (i.e.,  1948)  was  established  based  on  the  record-keeping 
practices  of  the  USAGE  for  locks  on  the  Mississippi  River.  This  stage  can  be 
characterized  as  a  stage  with  better  records  than  the  first  stage  that  are  needed  to 
compute  loading  cycles.  Generally,  the  available  information  in  this  stage  is  limited 
to  annual  lockages  and  annual  tonnage.  Therefore,  the  relationships  of  hardware 
cycles  and  lockage  cuts  as  functions  of  tonnage  that  were  developed  in  the  LPMS 
(third)  stage  can  be  used  to  estimate  hardware  cycles  and  lockages,  respectively. 
Alternatively,  the  number  of  hardware  cycles  as  a  function  of  time  and  the  number 
of  lockages  as  a  function  of  time  that  were  developed  in  the  LPMS  stage  can  be 
used  to  estimate  these  quantities  using  extrapolation. 


Loading  cycles  for  Stage  4  (present  to  planned  end  of  design  life  or 
rehabilitation) 

Stage  4  starts  from  the  present  and  ends  with  the  planned  end  of  design  (or 
rehabilitation)  life  of  a  lock.  The  models  in  this  stage  can  be  based  on  forecasting 
techniques.  The  scope  of  this  report  does  not  include  the  development  of 
forecasting  models.  However,  the  USAGE  General  Equilibrium  Model  (GEM) 
(USAGE  1994)  can  provide  forecasts  of  annual  tonnage  as  a  function  of  time 
based  on  a  set  of  input  variables.  The  predictions  of  traffic  volumes  are  expressed 
in  the  forms  of  low,  medium,  and  high  tonnage  predictions  as  functions  of  time. 
Therefore,  the  relationships  of  hardware  cycles  and  lockages  as  functions  of 
tonnage  that  were  developed  in  the  LPMS  (third)  stage  can  be  used  to  estimate 
hardware  cycles  and  lockage  cuts,  respectively,  for  this  stage. 


Probabilistic  Model 

The  objective  of  the  model  proposed  below  is  to  predict  the  number  of 
hardware  cycles  on  miter  gates  as  a  function  of  either  the  number  of  vessels 
passing  through  a  lock  or  the  number  of  lockage  cuts  that  occur  at  a  lock.  The 
model  accounts  for  the  following  factors:  (a)  navigation  traffic  volume,  (b)  traffic 
composition  in  terms  of  vessel  lengths,  (c)  length  and  capacity  of  navigation  locks, 
(d)  traffic  pattern  (in  terms  of  upstream/downstream  ratio),  and  (e)  non-vessel 
lockages  (loading  cycles  connected  with  passing  ice  and  debris). 
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Relationship  between  number  of  hardware  cycles  and  number  of 
vessels  (i.e.,  traffic  voiume) 

The  development  of  the  relationship  between  the  number  of  hardware  cycles 
and  the  number  of  vessels  (or  traffic  volume)  is  divided  into  the  following  four 
cases: 

a.  The  increase  in  the  number  of  hardware  cycles  due  to  cutting  of  long 
vessels,  assuming  the  traffic  is  in  one  direction. 

b.  The  effect  of  simultaneous  lockages  of  light  boats  or  recreational  vessels. 

c.  The  decrease  in  the  number  of  hardware  cycles  due  to  travel  in  opposite 
(upstream  and  downstream)  directions. 

d.  The  effect  of  environmental  conditions  (debris  and  ice  lockages). 

These  four  cases  are  discussed  below. 

Increase  in  number  of  hardware  cycles  due  to  cutting  of  long  vessels.  The 

maximum  length  (Imax)  of  a  vessel  that  can  be  locked  in  one  operation  of  a  lock  and 
the  cumulative  distribution  function  of  length  of  the  vessel  population  (Fl(1)) 
determine  the  number  of  needed  cuts.  The  length  of  the  vessel  population  is  a 
discrete  random  variable  with  the  following  cumulative  distribution  function: 


Fl(1)  =  2  Pjh)  for  i  =  1,2 . NL  (4) 

forlfil 


where  pdh)  =  probability  mass  value  for  a  vessel  of  a  length  /,  for  NL  possible 
discrete  vessel  lengths.  The  continuous  approximation  of  this  distribution  function 
is  an  integral  of  the  corresponding  probability  density  function 

I 

FJl)  =\f^{x)dx  (5) 

0 

The  arrival  of  vessels  in  one  direction  at  a  lock  can  be  assumed  to  follow  a 
Poisson  distribution  with  a  rate  2,  .  For  a  time  period  of  interest  (7),  which  is 
nonrandom,  the  mean  number  of  vessels  that  arrives  at  the  lock  during  the  time  T 
is 

N.  =  at  (6) 
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where  iVv  =  mean  number  of  vessels  arriving  in  time  T.  The  probability  (Pj)  that  a 
vessel  is  not  cut  in  two  or  more  parts  is  given  by 


Pi  —  Fldmax.) 


(7) 


The  probability  (P2)  that  a  vessel  is  cut  in  two  parts  is  given  by 

P2  =  FL^lmsn)  -  Fldmax)  (8) 

Analogously,  the  probability  (/**)  that  the  vessel  is  cut  in  ^  parts  is  given  by 


Pk  =  FL(kLj  -  Ft[(A-l)/n,»] 


(9) 


where  k  =  number  of  parts  into  which  vessels  can  be  cut.  In  general,  the 
probability  that  for  a  given  number  of  vessels  Ny  exactly  Nj  vessels  will  not  be  cut, 
N2  will  be  cut  in  two  parts,  and  so  on,  is  given  by  the  multinomial  distribution  with 
parameters  N^,  Pi,  P2,...,  Pk- 

The  probabilities  Pi,  P2,...,  Pk  should  satisfy  the  following  condition: 

k 

(10) 

i=l 

If  Nv  (the  number  of  vessels  that  arrive  at  a  lock  in  one  direction  during 
reference  time  period  7)  is  a  nonrandom  variable,  the  number  of  vessels  with  one 
or  more  parts  (i  =  parts)  due  to  cutting  has  the  following  mean  CNi),  and 

variance  {Var(Ni)),  respectively: 


Ni  =  N.Pi 


and 


Var(Ni)=N.Pi(l-Pi) 


(11) 

(12) 


The  number  of  vessels  which  are  cut  in  two  parts  has  the  following  mean  ( iVa )  and 
variance  (Var(N2)),  respectively: 


N2  =  N.P2 


(13) 
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and 


Var(N2)  =  N.P2(l-P2)  (14) 

The  number  of  vessels  which  are  cut  in  k  parts  has  the  following  mean  )  and 
variance  (Var(Nk)),  respectively: 


Nk  =  N.Pk  (15) 


and 

Var(Nic)  =  N.Pk(l-Pk)  (16) 

For  an  inbound  vessel  to  a  lock  with  its  gates  in  a  favorable  position  (i.e.,  a  fly 
entry  with  lock  chamber  pool  at  incoming  elevation)  that  needs  to  be  cut  into  / 
parts,  the  number  of  hardware  cycles  (iV^^c,  )  is  related  to  Ny  and  P,-  according  to 
Equation  3  as 


iV^^,=  (2i-l)N^P:  for/=l,2,...,^  (17) 

Therefore,  the  mean  of  the  total  number  of  hardware  cycles  ( Nfj^ )  can  be  related 
to  the  number  of  vessels  as 


Nhc  =  Ny{P2  +  3P2  +  5Ps+...+(2k-l)Pk)  (18) 


or 


Nhc  ~  Ny 


(SF;)+2P2  +  4Pj+...+2(^  — 1)P^ 
-  <=i 


(19) 


Using  Equation  10,  the  mean  total  hardware  cycles  can  be  expressed  as 


Nxic  ~  N  loc 


I  +  XO'-  i)Pi 


i=2 


(20) 
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The  variance  of  the  total  number  of  hardware  cycles  can  be  written  as 

Var(NHc)  =  - \fPil-P)  (21) 

i=l 


Decrease  in  number  of  hardware  cycles  due  to  simultaneous  lockages  of 
small  vessels.  It  is  common  in  operating  locks  to  simultaneously  service  a  number 
of  small  vessels,  as  is  the  case  for  light  boats  and  recreational  boats.  The  effect  of 
this  practice  on  hardware  cycles  is  a  reduction  in  its  total  number.  Therefore,  the 
above  approach  can  be  generalized  to  take  this  decrease  into  account.  Let  ps  be 
the  probability  that  a  given  boat  is  being  serviced  simultaneously  with  other  boats, 
then  the  mean  total  number  of  hardware  cycles  should  be  decreased  by  the  value 
ANhc 

(22) 

The  combination  of  Equations  19  and  22  produces 

Nhc  l+X2{i-])P,-2p,  (23) 

>■=2 

In  this  case,  the  following  condition  needs  to  hold 

i;pi+A=i  (24) 

i=I 

The  mean  of  the  total  number  of  hardware  cycles  according  to  Equation  23  is  a 
linear  function  of  the  number  of  vessels  for  trafSc  in  one  direction.  For  cases 
where  ^  the  increase  in  the  number  of  hardware  cycles  due  to  the  cutting 
of  long  vessels  prevails  over  the  decrease  in  the  number  of  hardware  cycles  due  to 
simultaneous  lockages  of  small  vessels,  and  vice  versa. 

In  the  LPMS  database,  the  number  of  lockage  cuts  (Nioc)  for  vessels  is  defined 
as 

(25a) 
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where  Ncuts  =  number  of  cuts,  then  Equation  25a  can  be  rewritten  using  Equations 
11,  13,  and  15  as 


(25b) 


where  Ncuts  takes  the  values  of  1,  2,  3, k,  in  which  1  corresponds  to  a  vessel 
without  a  cut.  In  this  case.  Equations  17,  22,  and  23  produce  the  following: 

(2«a) 

Equation  26a  can  be  rewritten  using  Equation  25b  to  produce  the  following 
relationship  between  mean  total  hardware  cycles  and  number  of  lockages  cuts: 

(2a>) 


Decrease  in  number  of  hardware  cycles  due  to  travel  in  opposite 
directions.  In  order  to  model  the  effect  of  two-direction  traffic  on  hardware 
cycles,  let  Xu  be  the  Poisson  rate  of  traffic  moving  upstream,  >id  be  the  Poisson 
arrival  rate  to  a  lock  for  traffic  moving  in  the  downstream  direction,  and  Nd  and 
Nvu  be  the  overall  number  of  vessels  arriving  to  a  lock  from  downstream  and 
upstream  directions  during  the  reference  period  T,  respectively.  If  dtd  is  the  service 
time  in  the  lock  for  a  vessel  in  the  downstream  traffic,  the  mean  decrease  in  the 
number  of  hardware  cycles  {ANncd)  for  a  steady-state  traffic  is 


ANnCd  —  NvdXup^  td 


(27) 


Since  N^d  +  N^  =  N,  the  downstream  number  of  vessels  Nd  =  Nv(J-a),  where  0  < 
ot  <  1,  Nvu  =  cNv,  and  a  =  fraction  of  traffic  moving  upstream.  Using  the  estimate 
Xu  =  Nvu/T,  Equation  27  can  be  expressed  as 


^Nncd 


a(l-a)Nldtd 

T 


(28) 


Analogously,  the  mean  decrease  in  the  number  of  hardware  cycles  due  to  the 
traffic  in  the  opposite  (upstream)  direction  is 


^NhCu 


a(l-a)Nldtu 

T 


(29) 
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where  Stu  is  the  service  time  in  the  lock  for  a  vessel  in  the  upstream  traffic,  and 
ANhcu  is  the  mean  decrease  in  the  number  of  hardware  cycles  in  the  upstream 
traffic.  The  sum  of  Equations  28  and  29  is  the  total  decrease  in  the  number  of 
hardware  cycles  {AN net ),  which  can  be  expressed  as 


Ai.r  _a(l-a)Nl(Stu  +  5u) 

^  JSHCt  — - 


(30) 


Equation  30  makes  sound  physical  sense.  For  example,  the  total  decrease  has  a 
maximum  value  when  the  upstream  traffic  is  equal  to  the  downstream  traffic  i.e., 
a  =  0.5.  Equation  30  can  be  re-written  as 


^Nnct 


a(l-a)Nl2a 

T 


(31) 


where  St  =  S  tup  +  St  down  /  2 .  Table  2  shows  estimates  of  hardware-cycle 

reduction  for  selected  values  of  a  and  N  using  a  =  0.5,  and  T  =  8,760  hours  (i.e., 
one  year). 


Table  2 
Estimates  of 
and  Ny  Using 

Hardware-Cycle  Reduction  for  Selected  Values  of  St 
a  -  0.5 

dt  (hours) 

For  Nv=^  ,000  vessels  per  year 

For  Ny  =  2,000  vessels  per 
year 

0.25 

14 

57 

0.5 

29 

114 

1 

57 

228 

2 

114 

457 

3 

171 

685 

4 

228 

913 

A  strong  seasonal  variation  in  the  number  of  vessels  can  result  in  considerable 
variation  of  a.  Thus,  in  this  case,  the  number  of  vessels  passing  through  the  lock 
cannot  be  considered  to  be  a  Poisson  process.  Under  this  restriction.  Equation  6 
needs  to  be  used  as  follows; 

a.  The  reference  time  period  in  this  case  can  be  taken  as  a  month,  i.e.,  T, 
where  i=  1,  2,...,  12. 

b.  The  fraction  of  upstream  traffic  is  a,,  which  is  related  to  the  /th  month. 
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c.  The  number  of  vessels  in  the  ith  month  is  Nvi,  so  that  if  Nv  is  the  annual 
number  of  vessels,  then 


(32) 

where  Pi  is  the  fraction  of  vessels  in  the  ith  month  from  the  traffic  of  a  year. 
Therefore,  the  following  condition  needs  to  be  satisfied: 


2;a='- 


(33) 


The  set  of  Pi  values  expresses  the  monthly  distribution  of  annual  traffic  in  both 
directions.  It  should  be  noted  that  neither  Pi  nor  a,  depends  on  the  absolute  value 
of  the  armual  number  of  vessels  TVy.  Thus,  in  this  case.  Equation  3 1  takes  on  the 
following  form: 


or 


2StNl f  X  2 
A  Nmt - - —  ^  (l  -  )Pi 


i=l 


^Nnct 


(34) 


(35) 


12 

where K='^a^{\ -  a^p\  is  a  coefficient  for  expressing  seasonal  variation  in 

i=I 

trafiBc  volume  and  direction. 


Hardware  cycles  due  to  non-vessel  lockages.  Some  of  the  hardware  cycles 
for  miter  gates  can  be  attributed  to  non-vessel  lockages  such  as  the  passing  of 
debris  and  ice  especialy  during  the  winter  months.  If  these  lockages  are  recorded 
in  the  LPMS  as  real  lockages  with  the  appropriate  numbers  of  hardware  cycles, 
they  can  be  taken  into  account  by  the  model  by  adding  a  term  to  Equation  23.  If 
these  lockages  are  not  recorded  in  the  LPMS,  then  their  hardware  cycles  can  be 
taken  into  account  by  adding  a  positive  constant  term  in  any  developed  regression 
models  to  fit  real  data  for  the  number  of  hardware  cycles  and  annual  tormage  as 
described  below. 
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Relationship  between  hardware  cycles  and  lockages  using 
regression  analysis 

In  this  section,  an  example  regression  model  is  developed  in  order  to  establish  a 
relationship  between  hardware  cycles  and  lockage  cuts  using  real  data,  e.g.,  based 
on  the  LPMS.  Equations  23  and  31  can  be  combined  to  obtain  the  following 
model  for  monthly  data  with  a  constant  term  {Constant)  that  corresponds  to  the 
hardware  cycles  associated  with  non- vessel  lockages  and  can  be  added  to  the 
model.  This  constant  would  not  be  needed  if  these  non- vessel  lockages  were 
recorded  in  the  LPMS  as  actual  lockages  with  the  appropriate  numbers  of 
hardware  cycles.  The  equation  for  this  model  is  given  by 


N  HC  ~ 


1=2 


a(l-a)N^2St 

f 


+  Constant 


(36a) 


The  mean  hardware  cycles  according  to  Equation  36a  can  be  expressed  in  terms  of 
Nioc  instead  of  Ny,  based  on  Equations  25b  and  26b  as 


=  *  Constant 


l'=l 


f  ^ 
V«=1 


(36b) 


or 


N„c^AN,or-BNl^  +  C 


(37) 


where  A,  B,  C  are  model  coefficients  with  B  >  0.  Figure  1  shows  a  scatter  diagram 
of  lockage  cuts  and  hardware  cycles  using  the  annual  data  that  are  provided  in 
Chapter  5.  The  scatter  diagram  shows  a  linear  relationship  between  hardware 
cycles  and  lockage  cuts,  i.e.,  B  =  0  in  Equation  37.  A  linear  regression  analysis 
was  then  performed.  For  the  case  B  =  0  and  C  =  0,  coefficient  A  has  the  same 
meaning  as  Kc  used  in  the  GEM  (USACE  1994),  i.e.,  the  mean  hardware  cycles 
per  lockage.  An  estimate  of  the  correlation  coefficient  was  obtained  to  be  0.999. 
The  estimates  of  the  coefficients  of  Equation  37  and  their  corresponding  standard 
errors  are  given  by 


A=  1.67898  ±0.0213 

(38) 

B  =  0  (not  significant) 

(39) 

-55.1587  ±  134.07  (not  significant) 

(40) 
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Figure  1 .  Annual  lockages  and  hardware  cycles 


Based  on  the  data  used  in  this  analysis,  only  the  coefficient  A  is  significant,  and 
should  be  kept  in  the  model.  The  developed  model  is  of  an  adequate  precision 
level  for  all  practical  purposes  since  the  standard  error  of  estimates  is  148.85  for  a 
sample  size  of  15  annual  values.  The  results  and  observations  provided  herein  are 
lock-specific.  For  data  obtained  from  other  locks,  with  different  patterns  of  traffic, 
the  significance  of  the  coefficients  A,  B,  and  C  can  be  different. 


T rend  analysis  of  annual  number  of  lockages 

The  objective  of  this  analysis  is  to  model  the  trend  of  the  annual  number  of 
lockages  for  selected  locks.  The  annual  number  of  lockages  is  considered  to 
follow  a  Poisson  distribution  with  the  following  mean 


Niccit)  =iVoexp(Ar  +  Br^  +  Ct^ +....) 


(41) 


where  t  is  the  time  in  years  counted  from  a  specified  year  (for  example  1948  for 
Locks  1 1  and  22,  and  1940  for  Lock  and  Dam  24).  Thus  t  =  0  for  1948  (or  1940), 
r  =  1  for  1949  (1941),  and  so  on.  The  coefficients  of  the  model ,  i.e..  No,  A,  B,  and 
C,  are  the  parameters  to  be  estimated  on  the  basis  of  curve  fitting  of  the  data.  In 
other  words,  the  sequence  of  annual  numbers  of  lockages  is  considered  to  be  a 
realization  of  a  nonhomogeneous  Poisson  process.  To  fit  the  model  of  Equation 
41  to  the  data,  a  regression  analysis  method  that  was  based  on  counts  with 
logarithmic  transformation  (Cox  and  Lewis  1966)  was  used.  The  significance  of 
the  coefficients  was  estimated  using  a  stepwise  regression  procedure.  Models  for 
Locks  24,  22,  and  1 1  were  fitted.  For  Lock  and  Dam  24,  the  following  estimates 
of  the  coefficients  and  their  standard  errors  were  obtained  using  the  data  in  Table  6 
for  the  years  1981  to  1993  and  the  data  provided  in  the  major  rehabilitation  report 


for  Lock  and  Dam  24  (USAGE  1993c): 

A(,=  1306.66  +  1.06  (42) 

A  =  0.050249  ±  0.00341 1  (in  years'^)  (43) 

B  =  0  (not  significant)  (44) 

C  =  -7.4844  X  lO'^  +  1.2260  x  lO"^  (in  years’^)  (45) 
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The  correlation  coefficient  for  the  model  is  0.955.  The  variance  of  the  residuals 
of  the  annual  lockages  is  80154.6.  The  corresponding  standard  error  for  the  model 
is  283.1.  The  accuracy  of  the  model  can  be  also  assessed  using  the  sum  of  fitted 
values  of  the  number  of  lockages  during  the  given  period.  The  sum  is  equal  to 
225132  which  does  not  differ  considerably  from  the  real  data  value  of  227147  (the 
difference  is  about  0.9%).  The  observed  and  fitted  (or  predicted  using  the  model) 
values  of  annual  numbers  of  lockages  are  given  in  Table  3.  The  results  are  also 
shown  in  Figure  2. 


Table  3 

Observed  and  Fitted  Annual  Number  of  Lockages  for 
Lock  and  Dam  24 


Year 

Observed 

Fitted 

Year 

Observed 

Fitted 

1940 

1718 

1307 

1971 

5250 

4964 

1941 

2212 

1374 

1972 

5584 

5104 

1942 

1409 

1445 

1973 

4914 

5242 

1943 

1050 

1519 

1974 

5157 

5375 

1944 

1119 

1597 

1975 

5058 

5503 

1945 

1213 

1678 

1976 

5386 

5625 

1946 

1817 

1763 

1977 

5235 

5741 

1947 

1705 

1853 

1978 

5668 

5849 

1948 

2178 

1946 

1979 

5778 

5949 

1949 

2551 

2043 

1980 

6601 

6040 

1950 

2258 

2144 

1981 

6788 

6121 

1951 

1990 

2249 

1982 

6319 

6193 

1952 

2027 

2357 

1983 

7108 

6253 

1953 

2335 

2470 

1984 

6595 

6302 

1954 

2435 

2587 

1985 

5099 

6338 

1955 

2981 

2707 

1986 

5262 

6363 

1956 

3146 

2831 

1987 

6523 

6373 

1957 

3288 

2959 

1988 

7131 

6370 

1958 

3386 

3090 

1989 

6787 

6354 

1959 

3698 

3224 

1990 

7743 

6324 

1960 

3408 

3362 

1991 

6963 

6280 

1961 

3365 

3502 

1992 

7228 

6222 

1962 

3801 

3644 

1993 

4189 

6150 

1963 

3827 

3789 

1994 

NA 

6064 

1964 

4084 

3935 

1995 

NA 

5965 

1965 

3585 

4083 

1996 

NA 

5854 

1966 

4078 

4231 

1997 

NA 

5729 

1967 

4211 

4379 

1998 

NA 

5594 

1968 

4688 

4527 

1999 

NA 

5447 

1969 

4344 

4675 

2000 

NA 

5290 

1970 

4938 

4820 

NA  =  not  applicable  or  not  available  | 
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Figure  2.  Observed  and  fitted  annual  number  of  lockages  for  Lock  and  Dam  24 


For  Lock  22  the  following  estimates  of  the  model  coefficients  and  their  standard 
errors  were  obtained  using  data  provided  by  the  USAGE: 


No  =  624.53  ±1.02 

(46) 

A  =  0.08392  ±  0.00256  (in  years'') 

(47) 

B  =  0  (not  significant) 

(48) 

G  =  -0.00002  ±  2.5x10-^  (in  years'^) 

(49) 

The  correlation  coefficient  for  this  model  is  0.996.  The  observed  annual  number 
of  lockages  and  the  fitted  ones  using  the  above  models  are  given  in  Table  4.  The 
results  are  also  shown  in  Figure  3. 

For  Lock  1 1  the  following  estimates  of  the  model  coefficients  and  their 
standard  errors  were  obtained  using  data  provided  by  the  USAGE: 


No  =  388.00  ±  1.04 

(50) 

A  =  0.10932  ±  0.00571  (in  years"') 

(51) 

B  =  -0.00147  ±  0.00017  (in  years'^) 

(52) 

C  =  0  (not  significant) 

(53) 

In  this  case,  the  correlation  coefficient  for  the  model  is  0.992.  The  observed  annual 
number  of  lockages  and  the  fitted  ones  using  the  above  models  are  given  in  Table  4. 
The  results  are  also  shown  in  Figure  3. 


Number  of  lockages  as  a  function  of  tonnage 

For  cases  where  data  on  the  annual  number  of  lockages  are  absent  but  the 
annual  tonnage  data  are  available,  the  relationship  between  annual  number  of 
lockages  and  tonnage  can  be  useful  to  estimate  the  number  of  lockages.  This 
relationship  can  be  obtained  on  the  basis  of  data  analysis  for  the  periods  where  both 
lockages  and  tonnage  information  is  available. 

The  model  development  requires  the  knowledge  of  the  number  of  lockages  and 
the  corresponding  tonnage  for  a  number  of  years.  For  each  record  value  of  annual 
number  of  lockages,  the  corresponding  annual  tonnage  value  is  needed.  This  model 
does  not  explicitly  account  for  recreational  boats  which  do  not  have  tonnage  values. 
Therefore,  the  4annual  number  of  lockages  Nioc  can  be  related  to  annual  tonnage  T„ 
(in  kilotons)  as 


Chapter  4  Methodology  for  Determining  Loading  Cyoles 


27 


(54) 


N^,  =  AT„+C 


where  the  coefficient  C  can  be  associated  with  recreational-boat  lockages,  or  with 
passing  ice  or  debris.  The  estimation  of  the  coefficients  for  this  model  using  Lock 
and  Dam  24  annual  data  for  the  period  1940  to  1993  produced  the  following 
estimates  for  the  coefficients:  C  =  1682.60  ±  92.13,  and  A  =  0.1432648  ± 
0.00424.  The  standard  error  of  estimates  for  the  model  is  395.96  based  on  54 
annual  values.  The  correlation  coefficient  between  annual  tonnage  and  number  of 
lockages  is  0.978. 


Table  4 

Observed  and  Fitted  Annual  Number  of  Lockages  for  Locks  22  and  1 1 

Year 

Lock  22 

1  Lock  11  1 

Observed 

Fitted 

Observed 

Fitted 

1948 

630 

625 

394 

388 

1949 

695 

680 

447 

432 

1950 

830 

738 

525 

480 

1951 

782 

803 

526 

531 

1952 

755 

872 

474 

587 

1953 

928 

948 

585 

646 

1954 

1015 

1029 

699 

709 

1955 

1159 

1116 

863 

776 

1956 

1197 

1210 

870 

847 

1957 

1298 

1309 

894 

921 

1958 

1466 

1417 

1092 

999 

1959 

1580 

1531 

1130 

1081 

1960 

1646 

1651 

1223 

1166 

1961 

1660 

1779 

1221 

1254 

1962 

1849 

1914 

1321 

1344 

1963 

2136 

2056 

1544 

1437 

1964 

2229 

2204 

1514 

1531 

1965 

2252 

2358 

1448 

1627 

1966 

2636 

2517 

1800 

1724 

1967 

2692 

2681 

1839 

1822 

1968 

2731 

2851 

1698 

1919 

1969 

3023 

3023 

1961 

2015 

1970 

3540 

3198 

2381 

2110 

1971 

3450 

3374 

2282 

2203 

1972 

3919 

3550 

2596 

2294 

1973 

3624 

3724 

2248 

2381 

1974 

3917 

3895 

2510 

2464 

1975 

3821 

4061 

2331 

2543 

1976 

3985 

4221 

2403 

2616 

1977 

4024 

4371 

2563 

2683 

1978 

4590 

4512 

2842 

2745 

1979 

4530 

4641 

2720 

2800 

1980 

5182 

4756 

3172 

2846 
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Estimating  the  cumuiative  number  of  hardware  cycies 

In  this  section,  the  cumulative  number  of  hardware  cycles  is  estimated  for  Lock 
and  Dam  24  over  the  period  1940  to  2000  for  demonstration  purposes.  The 
estimation  of  the  accumulated  number  of  hardware  cycles  can  be  based  on  (a) 
Equation  37  which  was  obtained  from  the  LPMS  data  for  the  period  1980  to  1992 
for  Lock  24,  and  (b)  the  estimated  annual  data  on  number  of  lockages  for  Lock  24 
for  the  time  period  from  the  start  of  life  of  1940  to  198 1  based  on  Equation  54. 

Using  Equation  41  for  the  annual  number  of  lockages  as  a  function  of  time  and 
Equation  34,  the  accumulated  number  of  hardware  cycles  were  predicted  to  the 
year  2000.  The  results  are  summarized  in  Table  5.  Table  5  includes  also  error 
estimates  based  on  the  actual  and  estimated  lockages  and  hardware  cycles.  These 
error  values  range  from  +0.4%  to  -2.9%. 


Table  5 

Total  of  Lockages  and  Hardware  Cycies  for  Lock  and  Dam  24 

Time  Period 

Actual  Lockages 
(Estimated  Lockages ;  %  Error) 

Actual  Hardware  Cycles 
(Estimated  Hardware 
Cycles :%  Error) 

1940  to  1979 

136,919 

(137,445  ;  +0.4%) 

NA 

(228,561,  NA) 

1980  to  1993 

90,279 

(87,686  ;  -2.9%) 

150,805 

(146,450 :  -2.9%) 

1994  to  2000 

NA 

(39,943  ;  NA) 

NA 

(66,677 :  NA) 

1940  to  2000 

NA 

(227,360  :  NA) 

NA 

(441,688  :  NA) 

1  NA  =  not  available 

It  should  be  noted  that  the  annual  number  of  lockages  can  also  be  predicted 
based  on  economic  forecasts  in  terms  of  annual  tonnage  values  using  for  example 
the  GEM  model  (USAGE  1994),  and  a  relationship  between  the  number  of 
lockages  and  tonnage  similar  to  Equation  54.  This  prediction  can  be  more 
accurate  than  using  Equation  41.  A  combination  of  the  two  approaches  can  also 
be  used. 
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5  Lock  and  Dam  24  -  Case 


Introduction 

The  objective  of  this  chapter  is  to  demonstrate  the  use  of  the  methodology  as 
described  in  Chapter  4  for  assessing  the  number  of  hardware  cycles  for  example 
miter  gates  at  a  navigation  lock.  The  miter  gates  of  Lock  and  Dam  24  were 
selected  for  this  purpose  and  are  the  sames  gates  used  in  Chapter  4  to  demonstrate 
some  aspects  of  the  methodology. 


Pool  and  Tail  water  Elevations 

The  hydraulic  records  for  Lock  and  Dam  24  were  obtained  for  the  period  1975 
to  1994.  The  daily  hydraulic  information  on  pool  and  tailwater  elevations  were 
obtained  from  the  USACE.  Figures  4a  and  4b  show  these  variations  on  a  daily 
basis  for  pool  and  tailwater  elevations,  respectively.  These  figures  demonstrate 
clearly  the  daily  variability  in  water  elevations,  hence  the  need  for  modeling  the 
problem  on  a  daily  basis,  not  monthly  nor  annually.  Figures  5a  and  5b  show  these 
variations  on  a  daily  basis  for  the  years  1975  to  1994  for  pool  and  tailwater 
elevations,  respectively. 

In  operating  a  lock  and  dam,  a  lockmaster  tries  to  maintain  a  pool  water  level 
by  adjusting  the  dam’s  gates.  However,  several  factors  contribute  to  the  variability 
in  the  pool  water  elevation;  e.g.,  increased  flow  levels  in  the  river  and  actions  by 
the  lockmasters  of  upstream  and  downstream  locks  to  control  pool  water  levels. 
Adjustments  to  a  dam’s  gates  result  in  changes  in  the  tailwater  elevation  as  well  as 
the  pool  water  elevation.  Therefore,  pool  and  tailwater  elevations  are  expected  to 
be  highly  correlated.  The  pool  water  elevation  for  Lock  and  Dam  24  was  plotted 
as  a  function  of  the  tailwater  elevation  in  Figure  6.  From  this  figure,  the  following 
three  regions  were  identified  and  are  described  in  subsequent  sections; 

a.  Low  tailwater  elevation,  <  439  ft. 

b.  Medium  tailwater  elevation,  439  ft  <H,<  445.062  ft. 

c.  High  tailwater  elevation,  Ht  >  445.062  ft. 
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Figure  4a.  Pool  water  elevations  in  1975  for  Lock  and  Dam  24 
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Figure  4b.  Tailwater  elevations  in  1993  for  Lock  and  Dam  24 


Day  of  the  Year 


Figure  5a.  Tailwater  elevations  for  Lock  and  Dam  24 
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Figure  6.  Water  elevations  for  1 975  to  1 994  -  Lock  and  Dam  24 


Low  tailwater  elevation 


It  is  observed  that  the  first  region  can  be  ch^acterized  with  a  high  level  of 
control  in  maintaining  the  pool  water  elevation  by  the  lockmaster.  The  pool  water 
elevation,  in  this  region,  is  maintained  at  a  constant  level  regardless  of  the  tailwater 
elevation.  Correlation  analyses  of  pool  and  tailwater  elevations  show  that  the 
correlation  level  is  not  significant.  Based  on  the  hydraulic  record,  statistical 
characteristics  of  the  pool  water  elevation  in  this  region  were  computed.  The 
mean  pool  water  elevation  is  448.89  ft,  and  the  standard  deviation  is  0.26  ft.  The 
sample  size  for  this  region  is  3579,  which  gives  a  standard  error  of  0.0044.  The 
maximum  and  minimum  pool  water  elevations  are  449.98  and  445.57  ft, 
respectively.  The  resulting  model  is  shown  in  Figure  7. 


Medium  tailwater  elevation 

The  region  of  medium  tailwater  elevation  shows  some  correlation  with  pool 
water  elevation  as  shown  in  Figure  6.  Therefore,  linear  regression  analysis  was 
performed  resulting  in  a  correlation  coefficient  of  -0.858.  The  following  linear 
prediction  model  of  pool  water  elevation  was  developed  for  this  region  for  Lock 
and  Dam  24: 

=  (-0^531396  ±  0.0068267)//,  +  (691.72094  +  3.01333)  (55) 

for  439  <//,<  445.062  ft) 

where  Hp  is  the  predicted  value  of  Hp.  Each  coefficient  in  the  model  is  provided 
with  its  +  standard  error.  The  standard  error  of  estimates  for  the  model  is 
0.561 1723  based  on  a  sample  size  in  this  region  of  2356.  Both  the  slope  and  the 
intercept  in  Equation  55  are  statistically  significant.  The  resulting  model  is  shown 
in  Figure  7. 


High  tailwater  elevation 

The  region  of  high  tailwater  elevation  shows  strong  correlation  with  pool  water 
elevation  as  shown  in  Figure  6.  Therefore,  linear  regression  analysis  was 
performed  resulting  in  a  correlation  coefficient  of  0.996.  The  following  linear 
prediction  model  of  pool  water  elevation  was  developed  for  this  region  for  Lock 
and  Dam  24: 

Hp  =  (1.001 18  ±  0.00001525)//,  for  H  >  445.062  ft  (56a) 
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Figure  7.  Model  for  water  elevations  for  1975  to  1994  -  Lock  and  Dam  24 


The  slope  coefficient  in  the  model  is  provided  with  its  ±  standard  error.  The 
intercept  for  this  model  was  determined  to  be  0.  The  standard  error  of  estimates 
for  the  model  is  0.23762  based  on  a  sample  size  in  this  region  of  1207.  The 
resulting  model  is  shown  in  Figure  7,  which  can  be  expressed  as 


H 


p 


448.89  +  0.0044 

(-0.5531396  ±  0.0068267)//, 
(691.72094  +  3.01333) 

(1.00118  +  0.00001525)//, 


+ 


//,  <439  ft 
439<//,  <445.062^3? 
//,  >  445.062^ 


(56b) 


The  design  pool  and  tailwater  elevations  are  449.0  ft  and  438.6  ft,  respectively. 
Equation  56b  can  be  normalized  with  respect  to  the  design  pool  and  tailwater 
elevations  to  obtain  the  normalized  pool  water  elevation  ( )  as  a  function  of  the 

normalized  tailwater  elevation  (/?„,)  as  follows: 


H. 


pn 


0.999755  +  0.0000097 
(-0.5403274 +  0.0066685)H,„  + 
(1.5405812  +  0.0067112) 

(0.97799  + 0.00001 48)//,„ 


<  1.000912 

1.000912  <  //„  <  1.0147332^ 
//„  >  1.0147332yr 


(56c) 


where  the  normalized  pool  and  tailwater  elevations  are  given,  respectively,  by 


and 


'  pn 


449.0 


H.- 


438.6 


(56d) 


(56e) 
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Hardware  Cycles 


The  computations  of  the  daily  hardware  cycles  were  based  on  the  data  obtained 
from  the  LPMS  and  use  of  the  method  described  in  Chapter  4.  The  daily  hardware 
cycles  were  computed  and  adjusted  for  ice  hardware  cycles.  The  adjusted  daily 
hardware  cycles  are  shown  in  Appendix  A.  The  adjustment  for  the  ice  lockages 
was  based  on  time-lapsed  videotapes  in  the  winter  months  of  1993-1994  for  Lock 
and  Dam  22  and  Lock  and  Dam  25  (Patev  1995).  The  videotapes  showed  63  and 
75  ice  lockages,  respectively,  over  periods  of  77  and  65  days,  respectively. 
Therefore,  one  ice  lockage  per  day  was  assumed  and  added  to  the  computed 
lockage  cuts  from  the  LPMS  for  the  months  of  January  and  February  of  each  year. 
Similarly,  one  ice  hardware  cycle  per  day  was  assumed  and  added  to  the  computed 
hardware  cycles  from  the  LPMS  for  the  months  of  January  and  February  of  each 
year. 

Using  the  definition  of  a  lockage  cut  as  the  process  of  passing  one  cut  of  a 
vessel  or  several  vessels  through  a  lock,  the  number  of  lockage  cuts  on  monthly 
and  yearly  bases  are  shown  in  Table  6.  The  number  of  hardware  cycles  on  a  daily 
basis,  and  on  monthly  and  yearly  bases  corrected  for  ice  lockages  are  shown  in 
Appendix  A  and  Table  7,  respectively.  The  number  of  lockage  cuts  are  shown  in 
Figures  8  and  9  on  monthly  and  yearly  bases.  The  number  of  hardware  cycles  are 
also  shown  in  Figures  10  and  1 1  on  monthly  and  yearly  bases.  The  number  of 
hardware  cycles  are  shown  in  Figure  12  on  a  daily  basis.  These  figures  indicate  the 
monthly,  seasonal,  and  yearly  variations  of  these  quantities.  The  intention  behind 
developing  these  plots  was  to  investigate  the  need  of  developing  hardware-cycle 
histograms  based  on  monthly,  seasonal,  or  yearly  parameters.  However,  as 
described  below,  the  computed  small  correlation  level  between  water  elevation  and 
number  of  hardware  cycles  allowed  the  aggregation  of  aU  water-elevation  records 
in  one  model  without  regard  to  monthly,  seasonal,  or  yearly  variations. 

Tailwater-Hardware  Cycles  Analysis 

The  objective  in  this  section  is  to  aggregate  the  number  of  hardware  cycles  that 
are  associated  with  the  same  tailwater  elevation  based  on  their  daily  records 
obtained  in  the  sections  above.  A  graphical  correlation  analysis  between  the 
hardware  cycles  and  tailwater  elevation  was  performed  as  shown  in  Figures  13  and 
14.  The  estimated  correlation  coefficient  is  0.155  which  is  small  indicating  that  all 
daily  records  can  be  aggregated  without  regard  to  monthly,  seasonal,  or  yearly 
variations.  The  number  of  hardware  cycles  that  are  associated  with  the  same 
tailwater  elevation  were  aggregated  to  obtain  a  histogram  as  shown  in  Figures  15 
and  16.  The  histogram  values  are  shown  in  Table  8.  Based  on  the  daily  records  of 
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tailwater  elevations  and  the  corresponding  hardware  cycles,  a  weighted  average  of 
water  elevation  was  computed  to  be  440.685  ft.  The  total  number  of  cycles  in  the 
entire  period  is  150,938.  The  standard  deviation'df  the  weighted  tailwater 
elevation  is  4.5817  ft.  The  maximum  and  minimum  tailwater  elevations  are  453.71 
and  433.34,  respectively. 


I  Table  6 

Summary  of  Lockages  from  1980  to  1994  for  Lock  and  Dam  24 


(Assumption;  Add  one  ice  cycle  per  day  during  Jan  and  Feb  for  all  years)  I 
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Table  7 

Summary  of  Hardware  Cycles  from  1980  to  1994  for  Lock  and  Dam  24 


1 1  (Assumption:  Add  one  ice  cycle  per  day  during  Jan  and  Feb  for  all  years) 
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Figure  8.  Summary  (3-dimensional)  of  lockages  for  1980  to  1994  -  Lock  and  Dam  24 
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Figure  10.  Summary  (3-dimensional)  of  hardware  cycles  for  1980  to  1994  -  Lock  and  Dam  24 


Figure  12.  Summary  (2-dimensional)  by  year  of  hardware  cycles  for  1980  to  1994  -  Lock  and  Dam  24 
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Figure  13.  Daily  tailwater  elevation  and  hardware  cycles  for  1985  and  1986  -  Lock  and  Dam  24 
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Figure  14.  Daily  tailwater  elevation  and  hardware  cycles  for  1980  to  1994  -  Lock  and  Dam  24 
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Table  8 

Data  for  Tailwater  Elevation  and  Hardware  Cycles  Hist 

ogram 

interval  Start  of 
Tailwater  (ft) 

Interval  End  of 
Tailwater  (ft) 

Mid-Interval 
of  Tailwater 
(ft) 

Count  of 
Hardware 
Cycles 

Fraction  of 
Hardware 
Cycles 
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0 

0 

TOTAL 

150938 

1 

The  resulting  histogram  can  be  used  to  select  a  probability  distribution  model. 
The  normal,  lognormal,  WeibuU,  and  Gumbel  probability  distributions  were 
considered  as  candidate  models.  The  chi-square  goodness-of-fit  test  (Ang  and 
Tang  1984)  was  used  to  select  the  best  model  among  the  candidate  distributions. 
Figure  17  shows  the  histogram  with  the  four  candidate  distributions.  By 
examining  Figure  17,  and  based  on'the  results  of  the  chi-square  test,  the  WeibuU 
distribution  can  be  considered  to  be  the  best  model  among  the  candidate  ones. 
The  WeibuU’s  cumulative  distribution  function  {Fh  (hi))  for  taUwater  elevation 
loading  is 
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Figure  1 7.  Histogram  of  tailwater  elevation  and  hardware-cycle  fraction  with  function  fits  to  data  for  Lock  and  Dam  24 


^h,(4)  =  1-exp 


a 

(57) 


where  Hu„in  =  minimum  tailwater  elevation  which  was  taken  as  433  ft,  Ht„:ax  = 
maximum  tailwater  elevation  which  was  taken  as  454  ft,  H,  =  tailwater  elevation 
loading  as  a  random  variable  with  a  given  value  of  ht,  and  a  (shape)  and  b  (scale) 
are  the  parameters  of  the  distribution.  The  shape  {a)  and  scale  Qj)  parameters 
were  estimated  using  the  best  linear  invariant  estimation  (Mann,  Shaffer,  and 
Singpurwalla  1974)  to  be  1.4909  and  0.3812,  respectively.  The  mean  and  variance 
based  on  these  parameters  are  440.23  ft  and  24.329  ft^,  respectively.  These 
moments  are  approximately  equal  to  the  moments  computed  as  the  weighted 
average  and  variance  of  tailwater  elevation  of  440.685  ft  and  20.99  ft^, 
respectively.  Equation  57  can  be  expressed  in  terms  of  the  normalized  water 
elevation 


Relationships  Among  Tonnage,  Lockages, 
Hardware  Cycles,  and  Time 

The  objective  of  this  section  is  to  study  the  relationships  for  Lock  and  Dam  24 
among  tonnage,  number  of  lockages,  number  of  hardware  cycles,  and  time. 

Figures  18, 19,  and  20  show  the  trend  of  tonnage,  lockages,  and  hardware  cycles. 
Figures  21, 22,  and  23  show  the  trend  of  the  ratio  of  tonnage  to  lockages,  the  ratio 
of  tonnage  to  hardware  cycles,  and  the  ratio  of  lockages  to  hardware  cycles, 
respectively.  The  last  set  of  figures  shows  the  scatter  diagrams  for  the 
relationships  between  lockages  and  tonnage  (Figure  24),  hardware  cycles  and 
tonnage  (Figure  25),  and  hardware  cycles  and  lockages  (Figure  26).  The 
relationship  between  hardware  cycles  and  lockages,  the  trend  relationship  of  the 
annual  number  of  lockages,  and  the  relationship  between  lockages  and  tonnage 
were  developed.  Then,  a  tonnage  forecast  model  was  developed  based  on  data 
obtained  using  the  GEM  (USAGE  1994)  for  Lock  and  Dam  24. 
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(Tonnage  tntons) 


Figure  22.  Trend  of  the  ratio  of  tonnage  to  hardware  cycles  for  Lock  and  Dam  24 


Year 


Figure  23.  Trend  of  the  ratio  of  lockages  to  hardware  cycles  for  Lock  and  Dam  24 
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Figure  24.  Tonnage  and  lockages  from  1940  to  1994  -  Lock  and  Dam  24 


Figure  25.  Tonnage  and  hardware  cycles  from  1980  to  1994  -  Lock  and  Dam  24 
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Figure  26.  Lockages  and  hardware  cycles  from  1980  to  1994  -  Lock  and  Dam  24 


Hardware  cycles  as  a  function  of  lockages 

According  to  Equation  36b,  a  regression  model  was  developed  for  the 
relationship  between  the  annual  hardware  cycles  and  the  annual  lockages  using  the 
data  from  the  LPMS  given  in  Tables  6  and  7.  According  to  Equations  37  to  40, 
the  model  can  be  expressed  as 

=  (1.678977  +  0.0213)  Ai,.  -  (0)iV^,  +  (-55.1587  + 134.07)  (58) 


Based  on  the  data  used  in  this  analysis,  only  the  coefficient  A  is  significant,  and 
should  be  kept  in  the  model.  However,  both  A  and  C  are  used  herein.  The 
developed  model  is  of  an  adequate  precision  level  for  all  practical  purposes,  since 
the  standard  error  of  estimates  is  148.85  for  a  sample  size  of  15  annual  values. 
The  results  and  observations  provided  herein  are  lock-specific.  For  data  obtained 
from  other  locks,  with  different  patterns  of  traffic,  the  significance  of  the 
coefficients  A,  B,  and  C  can  be  different.  The  model  of  Equation  58  and  the  data 
are  shown  in  Figure  1. 
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Trend  analysis  of  annual  number  of  lockages 

The  trend  of  the  annual  number  of  lockages  can  be  assumed  to  follow  a  Poisson 
distribution  with  a  time  variant  mean  as  given  by  Equation  41 .  The  coefficients  of 
the  model ,  i.e.,  No,  A,  B,  and  C,  need  to  be  estimated  on  the  basis  of  curve  fitting 
of  the  data  (Cox  and  Lewis  1966).  The  resulting  model  was  obtained  in  Equations 
42  to  45  to  be 

N!oc(t)  =(1306.66±1.06)exp((0.Q50249±0.003411)t+(-7.4844±1.2260)A:10“®r^)  (59) 

where  t  is  the  time  in  years  counted  from  a  specified  year  (for  example,  1940  for 
Lock  and  Dam  24).  Thus,  r  =  0  for  1940,  f  =  1  for  1941,  and  so  on.  The 
significance  of  the  coefficients  was  estimated  using  a  Stepwise  Regression 
procedure.  The  correlation  coefficient  for  the  model  is  0.955.  The  corresponding 
standard  error  is  283.1.  The  observed,  fitted  (or  predicted  using  the  model)  values 
of  annual  number  of  lockages  are  given  in  Table  3.  The  results  are  also  shown  in 
Figure  2.  The  details  for  the  development  of  this  model  are  given  in  Equations  42 
through  45. 


Number  of  lockages  as  a  function  of  tonnage 

For  cases  where  data  on  the  annual  number  of  lockages  are  absent  but  the 
annual  tonnage  data  are  available,  the  relationship  between  annual  number  of 
lockages  and  tonnage  can  be  useful  to  estimate  the  number  of  lockages  described 
in  Equation  54.  The  model  development  requires  the  values  of  annual  number  of 
lockages,  and  the  corresponding  annual  tonnage  values.  This  model  does  not 
explicitly  account  for  recreational  boats  which  do  not  have  tonnage  values. 
Therefore,  the  4annual  number  of  lockages  Nioc  can  be  related  to  annual  tonnage  T„ 
(in  kilotons)  for  Lock  and  Dam  24  for  the  period  1940  to  1993  as 

iVfc,  =  (0.1 432648  ±  0.00424)r„  +  (1 682.60  ±  92.1 3)  (60) 


where  the  constant  in  Equation  60  can  be  associated  with  recreational-boat 
lockages,  or  with  passing  ice  or  debris.  The  standard  error  of  estimates  for  this 
model  is  395.96  for  a  sample  of  size  54  annual  values.  The  correlation  coefficient 
between  annual  tonnage  and  number  of  lockages  is  0.978.  The  model  of  Equation 
60  and  the  data  are  shown  in  Figure  27. 

Tonnage  forecast  using  the  GEM 

The  GEM  was  used  to  obtain  tonnage  forecasts  for  Lock  and  Dam  24  for  the 
years  2000,  2010,  2020,  2030, 2040,  and  2050.  The  GEM  results  consist  of  low. 
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medium,  and  high  forecasts.  Table  9  shows  the  GEM  forecasts  for  Lock  and  Dam 
24.  The  relationship  of  annual  lockages  as  a  fiinction  of  the  annual  tonnage 
(Equation  60)  was  used  to  obtain  annual  lockages  for  Lock  and  Dam  24  as  shown 
in  Table  9.  Then  the  annual  hardware  cycles  as  a  function  of  lockages  (Equation 
58)  were  used  to  obtain  the  annual  hardware  cycles  for  Lock  and  Dam  24  as 
shown  in  Table  10.  The  forecasts  of  annual  tonnage,  lockages,  and  hardware 
cycles  are  shown  in  Figures  28,  29, 30,  respectively. 

The  GEM  tonnage  forecast  was  developed  in  1987  for  the  years  2000, 2010, 
2020, 2030,  2040,  and  2050.  Comparing  these  forecast  values  with  recently 
reported  “real”  values  as  given  in  Figure  18  shows  clearly  that  the  GEM  tonnage 
forecast  has  actually  underestimated  tonnage.  An  assessment  of  this  type  needs  to 
be  used  to  evaluate  and  possibly  revise  the  GEM. 
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Tables 

GEM  Forecasts  of  Tonnage  and  Computed  Lockages  for  Lock  and  Dam 
24 

Tonnage  (1000  tons) 

Number  of  Lockages  || 

Year 

High 

Medium 

Low 

High 

Medium 

Low 

2000 

41362 

41277 

41018 

7609 

7597 

7560 

2010 

41550 

41496 

41343 

7636 

7628 

7607 

2020 

41593 

41548 

41503 

7642 

7635 

7629 

2030 

41622 

41580 

41542 

7646 

7640 

7635 

2040 

41631 

41607 

41568 

7647 

7644 

7638 

2050 

41633 

41619 

41595 

7648 

7646 

7642 

Table  10  ~  — 

GEM  Forecasts  of  Tonnage  and  Computed  Hardware  Cycles  for  Lock 
and  Dam  24 


Tonnage  (1000  tons) 

Number  of  Hardware  Cycles 

Year 

High 

Medium 

Low 

High 

Medium 

Low 

2000 

41362 

41277 

41018 

12775 

12755 

12692 

2010 

41550 

41496 

41343 

12820 

12807 

12770 

2020 

41593 

41548 

41503 

12831 

12820 

12809 

2030 

41622 

41580 

41542 

12838 

12827 

12818 

2040 

41631 

41607 

41568 

12840 

12834 

12825 

2050 

41633 

41619 

41595 

12840 

12837 

12831 
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Annual  Lockages  ^  Annual  Tonnage  (1000  tons) 


41700 


Ycai 


5  28.  Forecast  of  annual  tonnage  for  Lock  and  Dam  24 


Figure  29.  Forecast  of  annual  lockages  for  Loc*  and  Dam  24 
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V 


Figure  30.  Forecast  of  annual  hardware  cycles  for  Lock  and  Dam  24 


Impact  of  Results  on  Fatigue  Reliability 
Assessment 

The  impact  of  the  results  of  this  study  on  fatigue  reliability  assessment  was 
investigated  in  two  forms,  (a)  its  impact  on  a  computed  water-head  differential 
(fid),  and  (b)  its  impact  on  a  computed  reliability  index  (P  ).  The  former  impact 
assessment  is  more  accurate  than  the  latter  one  due  to  approximations  employed  in 
the  latter  assessment.  The  results  reported  below  can  be  considered  as  a 
preliminary  assessment  of  the  effect  of  the  computed  water-elevation  and  hardware 
cycles  on  fatigue  reliability.  A  complete  assessment  of  this  impact  is  recommended 
for  a  future  study. 


Impact  on  water-head  differential 

The  daily  water  elevation  records  were  used  to  compute  water-head  differential 
(Hd),  for  which  a  histogram  was  developed  without  regard  to  load  cycles  as  shown 
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in  Figure  31.  The  mean  head  differential  is  8.267  ft;  the  standard  deviation  is 
4.7636  ft;  the  coefficient  of  variation  is  0.576;  and  the  standard  error  for  the  mean 
is  0.0565  ft.  By  accounting  for  the  hardware  cycles,  the  histogram  was  re¬ 
evaluated  based  on  Figure  16  and  Equation  56b  as  shown  in  Figure  32.  The 
weighted  mean  and  standard  deviation  of  the  water-head  differential  based  on  the 
hardware  cycles  are  7.624  ft  and  4.937  ft,  respectively.  The  coefficient  of 
variation  in  this  case  is  0.647.  A  USAGE  study  (USAEWES  1994)  provided 
estimates  of  mean  and  standard  deviation  for  water-head  differential  for  Lock  and 
Dam  24  of  9.600  ft,  and  4.167  ft,  respectively.  These  estimates  show  that  a  more 
rigorous  computation  of  the  water-head  differential  can  produce  a  lower  level  of 
head  differential  due  to  the  effect  of  hardware  cycles. 


Impact  on  reliability  index 

The  statistics  of  water-head  differential  were  used  in  the  USAGE  procedure  for 
computing  fatigue  reliability  of  a  vertical  beam  for  Lock  and  Dam  24  as  described 
in  USAGE  (1994).  The  results  are  shown  in  Table  1 1  that  demonstrate  the  effect 
of  water-head  differential  on  the  estimated  reliability  index.  The  results  are 
approximate  since  only  the  mean  and  standard  deviation,  not  the  complete 
probabilistic  characteristics  that  are  available,  of  the  water-head  differential  were 
used. 


Dsily  Hod  Diflcrtnti  Jl  (6) 


Figure  31.  Fraction  of  days  for  water-head  differential  from  1975  to  1994  -  Lock  and 
Dam  24 
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6  Recommendations 


As  a  result  of  using  the  LPMS,  the  development  of  the  probabilistic  model,  and 
performing  the  analysis,  several  recommendations  are  given  below. 

The  LPMS  entries  were  used  in  this  study  to  compute  hardware  cycles  that  are 
needed  for  fatigue-reliability  evaluation.  However,  the  computations  of  hardware 
cycles  require  a  significant  level  of  effort.  The  computation  of  hardware  cycles  can 
be  facilitated  by  adding  a  field  to  the  LPMS  to  keep  track  of  these  hardware 
cycles,  associate  times,  and  pool  and  tailwater  elevations.  Counting  hardware 
devices  with  timers  and  pressure  sensors  can  be  developed  to  perform  this 
function.  Alternatively,  current  fields  of  the  LPMS  can  be  improved  to  facilitate 
the  computations  of  hardware  cycles.  Also,  the  following  observations  are 
provided  as  suggestions  for  the  LPMS: 

a.  The  current  entry  and  exit  types  in  the  LPMS  do  not  necessarily  reflect  the 
turnback  type  if  it  was  delayed,  ie.,  not  immediate  to  an  entry  or  exit, 
respectively.  Depending  on  the  use  of  these  fields  in  their  current  forms, 
either  new  fields  should  be  developed  that  correct  for  the  delayed 
turnback  occurrence,  or  the  current  fields  should  be  revised. 

b.  The  fields  of  the  LPMS  need  to  be  logically  connected  in  order  to  prevent 
erroneous  entries. 

c.  Sometimes  several  vessel  records  were  entered  in  the  LPMS  as  separate 
lockages,  but  these  vessels  were  serviced  in  the  same  operation  of  opening 
and  closing  of  miter  gates.  The  LPMS  does  not  keep  track  of  these  cases, 
thereby  complicating  the  computation  of  hardware  cycles. 

d.  Ice  and  debris  lockages  are  not  included  in  the  LPMS.  The  practice  of 
record  keeping  needs  to  be  revised  to  require  the  inclusion  of  these 
lockages. 

e.  Other  operations  of  the  gates  for  service,  inspection,  or  performance 
evaluations  are  not  recorded  in  the  LPMS.  A  similar  action  to  item  d  is 
recommended  for  these  operations. 

In  general,  locks  and  dams  on  the  Mississippi  River  can  be  classified  into 
groups.  A  typical  lock  can  be  analyzed  from  each  group,  in  addition  to  analyzing 
several  locks  in  a  selected  group,  to  produce  a  complete  understanding  of 
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hardware  cycles  of  miter  gates.  Relationships  and  variability  among  the  groups 
and  within  a  group  can  then  be  studied  and  understood. 

The  GEM  tonnage  forecast  was  developed  in  1987  for  the  years  2000, 2010, 
2020, 2030, 2040,  and  2050.  Comparing  these  forecast  values  with  recently 
reported  “real”  values  shows  clearly  that  the  GEM  tonnage  forecast  has  actually 
underestimated  tonnage.  An  assessment  of  this  type  needs  to  be  used  to  evaluate 
and  possibly  revise  the  GEM.  Studies  in  this  area  are  needed  and  recommended. 

This  study  includes  a  preliminary  assessment  of  the  effect  of  the  computed 
water-elevation  and  hardware  cycles  on  fatigue  reliability.  A  complete  assessment 
of  this  impact  is  recommended  for  a  future  study. 
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CCi 

A 

a 
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Pi 
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CNC 

DR 

DY 

AMhc 
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a^ 
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EOLl 

E0L2 

ET 

M) 

Fl(1) 

GEM 

HC 

ht 

Ha 

Hp 

Hp„ 

Ht 

Htn 

Htmax 


Parameter  of  the  Weibull  distribution 
Monthly  fraction  of  upstream  traffic 
Regression  coefficient 
Fraction  of  traffic  moving  upstream 
Parameter  of  the  Weibull  distribution 
Regression  coefficient 
Reliability  index 

Fraction  of  vessels  in  a  month  from  the  traffic  of  a 
year 

Regression  coefficient 
Corrected  number  of  lockage  cuts 
Direction  of  lockage  (up  or  down) 

Day  of  shift 

Decrease  in  the  mean  of  total  hardware  cycles  due 
to  simultaneously  servicing  multiple  boats 
Average  service  time  in  the  lock  for  a  vessel 
Service  time  in  the  lock  for  a  vessel  in  doAvnstream 
traffic 

Service  time  in  the  lock  for  a  vessel  in  the  upstream 
traffic 

End  of  lockage  time  (24  hr)  first  cut 
End  of  lockage  time  (24  hr)  last  cut 
Entry  type 

Probability  density  function  of  vessel  length  L 
Cumulative  distribution  function  of  length  of  the 
vessel  population 
General  Equilibrium  Model 
Number  of  hardware  cycles 
Tailwater  elevation  value 
Water-head  differential 
Pool  elevation  (or  height)  of  water 
Normalized  pool  elevation  (or  height)  of  water 
Tailwater  elevation  (or  height) 

Normalized  tailwater  elevation  (or  height) 
Maximum  tailwater  elevation  (or  height) 
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Minimum  tailwater  elevation  (or  height) 

Predicted  value  of  Hp 
Number  of  cuts 

Coefficient  for  expressing  seasonal  variation  in 
traffic  volume  and  direction 
Mean  hardware  cycles  per  lockage 
Maximum  length  of  a  vessel  which  can  be  locked  in 
one  operation  of  a  lock 
Number  of  lockages 
Lock  number 

Lock  Performance  Monitoring  System 
Lockage  type 

Rate  of  vessel  arrival  at  a  lock 
Poisson  arrival  rate  to  a  lock  of  vessels  for 
downstream  traffic 

Poisson  arrival  rate  to  a  lock  of  vessels  for 
upstream  traffic 
Month  of  shift 
Model  coefficient 

Number  of  vessels  which  are  not  cut 
Number  of  vessels  which  are  cut  into  two  parts 
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Number  of  vessels  arriving  at  a  lock  in  time  T 
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Mean  number  of  vessels  arriving  at  a  lock  in  time  T 
from  downstream  direction 
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Mean  number  of  vessels  which  are  cut  into  two 
parts 

Mean  number  of  vessels  which  are  cut  into  k  parts; 
Mean  of  total  number  of  hardware  cycles 
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Mean  decrease  in  the  number  of  hardware  cycles  in 
the  downstream  traffic 


B2 


Appendix  B  Notation 


ANhCu 
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Pi 
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Ti 

TN 
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To 
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Var(N2) 

Var(NHc) 

Var(N^ 

VT 
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XT 

YR 


Mean  decrease  in  the  number  of  hardware  cycles  in 
the  upstream  traffic 

Mean  decrease  in  the  number  of  hardware  cycles 
for  two-direction  traffic 
Probability  mass  value  of  a  vessel  length  L 
Probability  that  a  given  boat  is  being  serviced 
simultaneously  with  other  boats 
Probability  that  a  vessel  is  not  being  cut 
Probability  that  a  vessel  is  being  cut  into  two  parts 
Probability  that  a  vessel  is  being  cut  into  k  parts 
Start  of  lockage  time  (24  hr)  1st  cut 
Time  in  years 
Reference  time  period 

Monthly  reference  time  period  for  i  =  1,  2, 12; 
Tonnage 
Annual  tonnage 
Mean  annual  tonnage 
Model  coefficient  for  mean  annual  tonnage 
Variance  of  number  of  vessels  which  are  not  cut 
Variance  of  number  of  vessels  which  are  cut  into 
two  parts 

Variance  of  total  number  of  hardware  cycles 
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Year  of  shift 


Appendix  B  Notation 


B3 


I  REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  this  collection  of  Information  is  estimated  to  average  1  hour  per  response.  Irwluding  the  time  for  reviewing  Instructions,  searching  existing  data  sources,  gathering  and  maintaining 
tfte  data  needed,  and  completing  and  reviewing  the  collection  of  information.  SerKi  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  mduding  suggestions 
for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  arKf  to  the 
Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20^. 

1 .  AGENCY  USE  ONLY  (Leave  blank) 

2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

December  1995  Final  report 

4.  TITLE  AND  SUBTITLE 

Loading  Cycle  for  the  Fatigue  Reliability  Analysis  of  Miter  Gates 

5.  FUNDING  NUMBERS 

6.  AUTHOR(S) 

Bilal  M.  Ayyub  Mark  P.  Kaminskiy 

Robert  C.  Patev  Mary  Ann  Leggett 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

BMA  Engineering,  Inc. 

14205  White  Water  Way,  Damestown,  MD  20878-3974; 

U.S.  Army  Engineer  Waterways  Experiment  Station 

3909  Halls  Ferry  Road,  Vicksburg,  MS  39180-6199 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

Technical  Report  ITL-95-12 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Corps  of  Engineers 

Washington,  DC  20314 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

111.  SUPPLEMENTARY  NOTES  | 

1  Available  from  National  Technical  Information  Service,  5285  Port  Royal  Road,  Springfield,  VA  22 1 6 1 . 

12a.  DISTRIBUTIONTAVAILABIUTY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

12b.  DISTRIBUTION  CODE 

13.  ABSTRACT  (Maximum  200  words) 


Miter  gates  at  navigation  locks  experience  loading  cycles  from  emptying  and  filling  of  a  lock’s  chamber  as  they  are 
opened  to  allow  traffic  through  the  locks.  Reliability  analysis  of  miter  gates  at  navigation  locks  requires  definition  of 
(a)  nonperformance  modes,  (b)  loads,  (c)  structural  strength,  and  (d)  methods  of  reliability  analysis.  Due  to  the  cyclic  loading 
nature  of  miter  gates,  the  fatigue  of  critical  details  requires  examination  using  reliability  methods.  The  assessment  of  fatigue 
reliability  of  these  details  as  a  function  of  time  requires  the  knowledge  of  strength,  stress  ranges,  and  loading  cycles  for  these 
details. 

Prediction  of  loading  cycles  on  miter  gates  for  use  in  the  assessment  of  fatigue  reliability  for  miter  gates  is  described. 
Correlation  of  field  data  for  pool  and  tailwater  elevations  and  barge  traffic  to  form  a  loading  histogram  to  be  utilized  to  better 
predict  the  loading  history  of  miter  gates  is  explained. 


1  1 

|l4.  SUBJECT  TERMS 

Reliability  Lockages  Tailwater 

Miter  Gates  Locks  Pool  water 

Fatigue  Hardware  Cycles  Loading  histogram 

1  Rehability  Index 

15.  NUMBER  OF  PAGES 

91 

16.  PRICE  CODE 

17.  SECURITY  CLASSIRCATION 
OF  REPORT 

1  UNCLASSIFIED 

18.  SECURITY  CLASSIRCATION 
OF  THIS  PAGE 

UNCLASSIFIED 

19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

NSN  7540.01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  Z39>18 
298-102 


